Frontiers in plant RNA research in ICAR2023: from lab to innovative agriculture.

The recent growth in global warming, soil contamination, and climate instability have widely disturbed ecosystems, and will have a significant negative impact on the growth of plants that produce grains, fruits and woody biomass. To conquer this difficult situation, we need to understand the molecular bias of plant environmental responses and promote development of new technologies for sustainable maintenance of crop production. Accumulated molecular biological data have highlighted the importance of RNA-based mechanisms for plant stress responses. Here, we report the most advanced plant RNA research presented in the 33rd International Conference on Arabidopsis Research (ICAR2023), held as a hybrid event on June 5-9, 2023 in Chiba, Japan, and focused on "Arabidopsis for Sustainable Development Goals". Six workshops/concurrent sessions in ICAR2023 targeted plant RNA biology, and many RNA-related topics could be found in other sessions. In this meeting report, we focus on the workshops/concurrent sessions targeting RNA biology, to share what is happening now at the forefront of plant RNA research.

Fluorescent In Situ Detection of Small RNAs in Plants Using sRNA-FISH.

Plant small RNAs are 21-24 nucleotide, noncoding RNAs that function as regulators in plant growth and development. Colorimetric detection of plant small RNAs was made possible with the introduction of locked nucleic acid probes. However, fluorescent detection of plant small RNAs has been challenging due to the high autofluorescence from plant tissue. Here we report a fluorescent in situ detection method for plant small RNAs. This method can be applied to most plant samples and tissue types and also can be adapted for single-molecule detection of small RNAs with super-resolution microscopy.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation.

KEY MESSAGE: Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

Plant non-coding RNAs: The new frontier for the regulation of plant development and adaptation to stress.

Most plant transcriptomes constitute functional non-coding RNAs (ncRNAs) that lack the ability to encode proteins. In recent years, more research has demonstrated that ncRNAs play important regulatory roles in almost all plant biological processes by modulating gene expression. Thus, it is important to study the biogenesis and function of ncRNAs, particularly in plant growth and development and stress tolerance. In this review, we systematically explore the process of formation and regulatory mechanisms of ncRNAs, particularly those of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, we provide a comprehensive overview of the recent advancements in ncRNAs research, including their regulation of plant growth and development (seed germination, root growth, leaf morphogenesis, floral development, and fruit and seed development) and responses to abiotic and biotic stress (drought, heat, cold, salinity, pathogens and insects). We also discuss research challenges and provide recommendations to advance the understanding of the roles of ncRNAs in agronomic applications.

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants.

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

Harnessing the potential of non-coding RNA: An insight into its mechanism and interaction in plant biotic stress.

Plants have developed diverse physical and chemical defence mechanisms to ensure their continued growth and well-being in challenging environments. Plants also have evolved intricate molecular mechanisms to regulate their responses to biotic stress. Non-coding RNA (ncRNA) plays a crucial role in this process that affects the expression or suppression of target transcripts. While there have been numerous reviews on the role of molecules in plant biotic stress, few of them specifically focus on how plant ncRNAs enhance resistance through various mechanisms against different pathogens. In this context, we explored the role of ncRNA in exhibiting responses to biotic stress endogenously as well as cross-kingdom regulation of transcript expression. Furthermore, we address the interplay between ncRNAs, which can act as suppressors, precursors, or regulators of other ncRNAs. We also delve into the regulation of ncRNAs in response to attacks from different organisms, such as bacteria, viruses, fungi, nematodes, oomycetes, and insects. Interestingly, we observed that diverse microorganisms interact with distinct ncRNAs. This intricacy leads us to conclude that each ncRNA serves a specific function in response to individual biotic stimuli. This deeper understanding of the molecular mechanisms involving ncRNAs in response to biotic stresses enhances our knowledge and provides valuable insights for future research in the field of ncRNA, ultimately leading to improvements in plant traits.

Nonsense-mediated mRNA decay pathway in plants under stress: general gene regulatory mechanism and advances.

MAIN CONCLUSION: Nonsense-mediated mRNA decay in eukaryotes is vital to cellular homeostasis. Further knowledge of its putative role in plant RNA metabolism under stress is pivotal to developing fitness-optimizing strategies. Nonsense-mediated mRNA decay (NMD), part of the mRNA surveillance pathway, is an evolutionarily conserved form of gene regulation in all living organisms. Degradation of mRNA-bearing premature termination codons and regulation of physiological RNA levels highlight NMD's role in shaping the cellular transcriptome. Initially regarded as purely a tool for cellular RNA quality control, NMD is now considered to mediate various aspects of plant developmental processes and responses to environmental changes. Here we offer a basic understanding of NMD in eukaryotes by explaining the concept of premature termination codon recognition and NMD complex formation. We also provide a detailed overview of the NMD mechanism and its role in gene regulation. The potential role of effectors, including ABCE1, in ribosome recycling during the translation process is also explained. Recent reports of alternatively spliced variants of corresponding genes targeted by NMD in Arabidopsis thaliana are provided in tabular format. Detailed figures are also provided to clarify the NMD concept in plants. In particular, accumulating evidence shows that NMD can serve as a novel alternative strategy for genetic manipulation and can help design RNA-based therapies to combat stress in plants. A key point of emphasis is its function as a gene regulatory mechanism as well as its dynamic regulation by environmental and developmental factors. Overall, a detailed molecular understanding of the NMD mechanism can lead to further diverse applications, such as improving cellular homeostasis in living organisms.

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

The plant siRNA landscape.

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

PRMD: an integrated database for plant RNA modifications.

The scope and function of RNA modifications in model plant systems have been extensively studied, resulting in the identification of an increasing number of novel RNA modifications in recent years. Researchers have gradually revealed that RNA modifications, especially N6-methyladenosine (m6A), which is one of the most abundant and commonly studied RNA modifications in plants, have important roles in physiological and pathological processes. These modifications alter the structure of RNA, which affects its molecular complementarity and binding to specific proteins, thereby resulting in various of physiological effects. The increasing interest in plant RNA modifications has necessitated research into RNA modifications and associated datasets. However, there is a lack of a convenient and integrated database with comprehensive annotations and intuitive visualization of plant RNA modifications. Here, we developed the Plant RNA Modification Database (PRMD; http://bioinformatics.sc.cn/PRMD and http://rnainformatics.org.cn/PRMD) to facilitate RNA modification research. This database contains information regarding 20 plant species and provides an intuitive interface for displaying information. Moreover, PRMD offers multiple tools, including RMlevelDiff, RMplantVar, RNAmodNet and Blast (for functional analyses), and mRNAbrowse, RNAlollipop, JBrowse and Integrative Genomics Viewer (for displaying data). Furthermore, PRMD is freely available, making it useful for the rapid development and promotion of research on plant RNA modifications.

RIP-Seq analysis of non-PPR chloroplast editing factors reveals broad RNA interactions and enrichment of less efficiently translated RNAs by OZ1 and ORRM1 complexes.

C-to-U RNA editing in angiosperm chloroplasts requires a large suite of proteins bound together in the editosome. The editosome is comprised of PPR proteins, RIP/MORFs, OZ proteins, and ORRM proteins that physically interact in high molecular weight complexes. The specific functions of non-PPR editing factors in the editosome are unclear, however, specific subsets of editing sites are affected by absence of non-PPR editing factors. Unlike the PPR components of editosomes that have predictable nucleotide specificities, domains present in non-PPR editing factors make RNA associations difficult to predict. In this study, chloroplast extracts were isolated from juvenile maize seedlings. RNAs were immunoprecipitated using polyclonal antibodies targeting non-PPR editing factors RIP9, OZ1, and ORRM1. RNA libraries from duplicate experiments were compared. RIP9 was associated with most of the non-ribosomal RNA content of chloroplasts, consistent with a general binding function to PPR L-motifs and tethering of large ribonucleoprotein complexes. The breadth of RNA associations was greater than predicted and include mRNAs without predicted editing sites, tRNA sequences, and introns. OZ1 and ORRM1 were associated with a highly similar pool of RNAs that have a bias toward lower translational efficiency values in mature chloroplasts. Lower translational efficiency was also associated with the pool of edited RNAs compared to RNAs without editing sites. The unexpected breadth of interactions by non-PPR editing factors suggests the editosome is large, diverse, and associated with RNAs with lower relative translational efficiency in mature chloroplasts.

Profiling of RNA editing events in plant organellar transcriptomes with high-throughput sequencing.

RNA editing is a crucial post-transcriptional modification process in plant organellar RNA metabolism. rRNA removal-based total RNA-seq is one of the most common methods to study this event. However, the lack of commercial kits to remove rRNAs limits the usage of this method, especially for non-model plant species. DSN-seq is a transcriptome sequencing method utilizing duplex-specific nuclease (DSN) to degrade highly abundant cDNA species especially those from rRNAs while keeping the robustness of transcript levels of the majority of other mRNAs, and has not been applied to study RNA editing in plants before. In this study, we evaluated the capability of DSN-seq to reduce rRNA content and profile organellar RNA editing events in plants, as well we used commercial Ribo-off-seq and standard mRNA-seq as comparisons. Our results demonstrated that DSN-seq efficiently reduced rRNA content and enriched organellar transcriptomes in rice. With high sensitivity to RNA editing events, DSN-seq and Ribo-off-seq provided a more complete and accurate RNA editing profile of rice, which was further validated by Sanger sequencing. Furthermore, DSN-seq also demonstrated efficient organellar transcriptome enrichment and high sensitivity for profiling RNA editing events in Arabidopsis thaliana. Our study highlights the capability of rRNA removal-based total RNA-seq for profiling RNA editing events in plant organellar transcriptomes and also suggests DSN-seq as a widely accessible RNA editing profiling method for various plant species.

A universal protocol for high-quality DNA and RNA isolation from diverse plant species.

Next-generation sequencing demands high-quality nucleic acid, yet isolating DNA and RNA is often challenging, particularly from plant tissues. Despite advances in developing various kits and reagents, these products are tailored to isolation of nucleic acid from model plant tissues. Here we introduce a universal lysis buffer to separate nucleic acid from various plant species, including recalcitrant plants, to facilitate molecular analyses, such as quantitative PCR (qPCR), transcriptomics, and whole-genome sequencing (WGS). The protocol is a modification of the original CTAB methods, which leads to nucleic acid isolation from many plant species, including monocots and eudicots. The lysis buffer consists of hexadecyltrimethylammonium bromide (CTAB), sodium chloride (NaCl), Tris base, ethylenediaminetetraacetic acid (EDTA) and ß-mercaptoethanol (ßME). The modified CTAB method enables the isolation of nucleic acid from small amounts of plant tissues (e.g., 15-100 mg) in a timely manner, which is well-suited for a large number of samples and also when adequate sample collection is a limiting factor. The protocol isolates not only DNA from various plant species but also RNA. This makes it highly effective for molecular analyses compared to previously described CTAB methods optimised for DNA isolation. The appropriate concentration of the components enables high-quality DNA and RNA isolation from plant tissues simultaneously. Additionally, this protocol is compatible with commercially available columns. For DNA and RNA to be qualified for next-generation sequencing platforms, the protocol is supplemented with columns to purify either DNA or RNA from the same tissue to meet high standards for sequencing analyses. This protocol provides an ideal approach to overcome potential obstacles in isolating high-quality DNA or RNA from a wide range of plant species for downstream molecular analysis.

PhasiHunter: a robust phased siRNA regulatory cascade mining tool based on multiple reference sequences.

SUMMARY: In recent years, phased small interfering RNA has been found to play crucial roles in many biological processes in plants. However, efficiently predicting phasiRNA regulatory cascades with computational methods is still challenging. Here, we introduce PhasiHunter, a phasiRNA regulatory network prediction tool that has several distinctive features compared to existing tools: (i) PhasiHunter employs two major phasiRNA prediction algorithms, namely phase score and hypergeometric distribution-based methods, to ensure the integrity and accuracy of prediction; (ii) PhasiHunter can identify phasiRNAs and their regulatory networks based on multiple reference sequences and the predicted results can be automatically integrated; (iii) PhasiHunter can efficiently identify the phasiRNAs generated through alternative splicing events; and (iv) the excellent data structure and parallel computing architecture allow PhasiHunter to predict phasiRNAs and their regulatory pathways with high efficiency. AVAILABILITY AND IMPLEMENTATION: PhasiHunter is an open-source tool that is available at https://github.com/HuangLab-CBI/PhasiHunter.

Small RNA sequencing provides insights into molecular mechanism of flower development in Rhododendron pulchrum Sweet.

Rhododendron pulchrum sweet, a member of the Ericaceae family possessing valuable horticultural properties, is widely distributed in the temperate regions. Though serving as bioindicator of metal pollution, the molecular mechanism regulating flowering in R. pulchrum is very limited. Illumina sequencing was performed to identify critical miRNAs in the synthesis of flavonoids at different developmental stages. Totally, 722 miRNAs belonging to 104 families were screened, and 84 novel mature miRNA sequences were predicted. The miR166, miR156, and miR167-1 families were dominant. In particular, 126 miRNAs were significantly differentially expressed among four different flowering stages. Totally, 593 genes were differentially regulated by miRNAs during the flower development process, which were mostly involved in "metabolic pathways", "plant hormone signal transduction", and "mitosis and regulation of biosynthetic processes". In pigment biosynthesis and signal transduction processes, gra-miR750 significantly regulated the expression of flavonoid 3',5'-hydroxylase; aof-miR171a, aof-miR171b, aof-miR171c, cas-miR171a-3p, and cas-miR171c-3p could regulate the expression of DELLA protein; aof-miR390, aof-miR396b, ath-miR3932b-5p, cas-miR171a-3p, aof-miR171a, and aof-miR171b regulated BAK1 expression. This research showed great potentials for genetic improvement of flower color traits for R. pulchrum and other Rhododendron species.

The plant noncoding transcriptome: a versatile environmental sensor.

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Plant immunity suppressor SKRP encodes a novel RNA-binding protein that targets exon 3' end of unspliced RNA.

The regulatory roles of RNA splicing in plant immunity are emerging but still largely obscure. We reported previously that Phytophthora pathogen effector Avr3c targets a soybean protein SKRP (serine/lysine/arginine-rich protein) to impair soybean basal immunity by regulating host pre-mRNA alternative splicing, while the biochemical nature of SKRP remains unknown. Here, by using Arabidopsis as a model, we studied the mechanism of SKRP in regulating pre-mRNA splicing and plant immunity. AtSKRP confers impaired plant immunity against Phytophthora capsici and associates with spliceosome component PRP8 and splicing factor SR45, which positively and negatively regulate plant immunity, respectively. Enhanced crosslinking and immunoprecipitation followed by high-throughput sequencing (eCLIP-seq) showed AtSKRP is a novel RNA-binding protein that targets exon 3' end of unspliced RNA. Such position-specific binding of SKRP is associated with its activity in suppressing intron retention, including at positive immune regulatory genes UBP25 and RAR1. In addition, we found AtSKRP self-interact and forms oligomer, and these properties are associated with its function in plant immunity. Overall, our findings reveal that the immune repressor SKRP is a spliceosome-associated protein that targets exon 3' end to regulate pre-mRNA splicing in Arabidopsis.

Identification of microRNAs involved in ectomycorrhizal formation in Populus tomentosa.

The majority of woody plants are able to form ectomycorrhizal (ECM) symbioses with fungi. During symbiotic development, plants undergo a complex re-programming process involving a series of physiological and morphological changes. MicroRNAs (miRNAs) are important components of the regulatory network underlying symbiotic development. To elucidate the mechanisms of miRNAs and miRNA-mediated mRNA cleavage during symbiotic development, we conducted high-throughput sequencing of small RNAs and degradome tags from roots of Populus tomentosa inoculated with Cenococcum geophilum. This process led to the annotation of 51 differentially expressed miRNAs between non-mycorrhizal and mycorrhizal roots of P. tomentosa, including 13 novel miRNAs. Increased or decreased accumulation of several novel and conserved miRNAs in ECM roots, including miR162, miR164, miR319, miR396, miR397, miR398, novel-miR44 and novel-miR47, suggests essential roles for these miRNAs in ECM formation. The degradome analysis identified root transcripts as miRNA-mediated mRNA cleavage targets, which was confirmed using real-time quantitative PCR. Several of the identified miRNAs and corresponding targets are involved in arbuscular mycorrhizal symbioses. In summary, increased or decreased accumulation of specific miRNAs and miRNA-mediated cleavage of symbiosis-related genes indicate that miRNAs play important roles in the regulatory network underlying symbiotic development.

An mTRAN-mRNA interaction mediates mitochondrial translation initiation in plants.

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that "mitochondrial translation factors" mTRAN1 and mTRAN2 are land plant-specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)-like RNA binding proteins of the mitoribosomal "small" subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5' regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

RNA pseudouridine modification in plants.

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.